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Abstract 

In the framework of pionless nucleon-nucleon effective field theory we study different approxi- 
mation schemes for the nuclear many body problem. We consider, in particular, ladder diagrams 
constructed from particle-particle, hole-hole, and particle-hole pairs. We focus on the problem of 
finding a suitable starting point for perturbative calculations near the unitary limit (kpa) — > oo 
and (kpr) — > 0, where kp is the Fermi momentum, a is the scattering length and r is the effective 
range. We try to clarify the relationship between different classes of diagrams and the large g and 
large D approximations, where g is the fermion degeneracy and D is the number of space-time 
dimensions. In the large D limit we find that the energy per particle in the strongly interacting 
system is 1/2 the result for free fermions. 
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I. INTRODUCTION 

The nuclear many-body problem is of fundamental importance to nuclear physics [lj. The 
traditional approach to the many-body problem is based on the assumption that nucleons 
can be treated as non-relativistic point-particles interacting mainly via two-body potentials. 
Three-body potentials, relativistic effects, and non-nucleonic degrees of freedom are assumed 
to give small corrections. The many-body Schrodinger equation is solved using a variety of 
methods that involve both variational and numerical aspects. 

Over the last several years an alternative approach to nuclear physics based on effective 
field theory (EFT) methods has been applied successfully to the two and three-nucleon sys- 
tems EFT methods have the advantage of being directly connected to QCD, 
and of providing a framework which is amendable to systematic improvements. The appli- 
cation of EFT to nuclear systems is complicated by the appearance of anomalously small 
energy scales. In the two-body system these small scales are reflected by the large neutron- 
neutron and neutron-proton scattering lengths, a rm ( 1 S'o) — —18 fm and a np ( 3 5'i) ~ 5 fm, 
respectively. Effective field theories capable of describing systems with anomalously large 
scattering lengths require summing an infinite number of Feynman diagrams at leading order 

a. 

In this work we wish to study the EFT approach to the nuclear many body problem 
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121 ]. We will focus on the equation of state of pure neutron matter at low 



to moderate density, a problem that is of relevance to the structure of neutron stars. The 
neutron matter problem has the theoretical advantage that there are no three-body forces 
at leading order. As in the two-body system the main obstacle is the large scattering length. 
If the scattering length was small the equation of state and other quantities of interest 
could be expanded in (kpct), where hp is the Fermi momentum. This is the standard low 
density expansion for a hard sphere Fermi gas which was studied by Huang, Lee and Yang 
in the 1950 's and rederived in the EFT context by Hammer and Furnstahl 0H,Q- In 
real nuclear matter, however, \kpa\ ^> 1 and the perturbative low density expansion is not 
useful. 

An interesting system that illustrates the difficulties of the nuclear matter problem is 
a dilute liquid of non-relativistic spin 1/2 fermions interacting via a short range potential 
with infinite scattering length. In this case the parameters that characterize the many body 
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problem are either infinite or zero, \kpa\ — > oo and \kpr\ — > 0, where a,r are the scattering 
length and the effective range. Dimensional analysis implies that the equation of state is of 
the form 

A ? 52m' 1 J 

where £ is dimensionless number. For free fermions £ = 1, but for strongly correlated 
fermions the theory contains no obvious expansion parameter and the determination of £ 
is a difficult non-perturbative problem. Recent interest in this problem has been fueled by 
experimental advances in creating cold, dilute gases of fermionic atoms tuned to be near a 
Feshbach resonance These experiments are beginning to yield results for the equation 
of state of non-relativistic fermions in the limit \kpa\ — > oo. 

A plausible strategy for investigating neutron matter is to start from a numerical or 
variational solution of the "unitary limit" system and to include corrections due to the finite 
effective range, explicit pion degrees of freedom, many body forces, etc. perturbatively. 
Recent numerical studies of many body systems with a large scattering length can be found 
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In this work we study analytic many body approximations that could 
be used as the starting point for a theory of neutron matter based on EFT interactions. 
We focus on ladder diagrams built from particle-particle or particle-hole bubbles and study 
whether these approximations can be consistently renormalized and yield a stable \kpa\ — > oo 
limit. We also examine the relationship of these approaches to the large g and large D limits, 
where g is the number of fermion fields and D is the number of space-time dimensions. 

II. PARTICLE-PARTICLE LADDER DIAGRAMS 

We will consider non-relativistic fermions governed by an effective lagrangian of the form 

(2) 



C eS = iVt \^id + ^JN- y (iV f iV) 2 + [(jVJV)t(jVV N) + h.c. 

where V=V — V- We have not displayed terms with higher derivatives or more powers 
of the fermion field, including two-derivative terms that act in the p-wave channel. The 
parameters Cq and C2 are related to the s-wave scattering length and the effective range. 
In the power divergence subtraction (PDS) scheme the relationship is given by 2^1 

£5, = -^, C^^-MM^ (3) 
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FIG. 1: Particle-particle ladder diagrams (left panel) and particle-hole ring diagrams (right panel) 
in the effective field theory. 



where p is the renormalization scale. The advantage of the PDS scheme is that the Cj are of 
natural size even if the scattering length is large. Dimensional regularization with minimal 
subtraction gives Cq = (4na)/M which is unnatural in the limit a — > oo. 

We are interested in the energy density and pressure of a many body system with baryon 
density p. At zero temperature the density is related to the Fermi momentum via p = 
gkp/(6ir 2 ), where g = 2 is the degeneracy factor. The free propagator is given by 

G (kU - <W [ ko _ k2/2M + ie + ko _ k 2 /2M _ ie ) . ( 4 ) 

and describes two types of excitations, holes with momentum k < k F and particles with 
k > kp. Using this propagator and the vertices from equ. (J2J) we can compute the en ergy 
per particle as a perturbative expansion in (k F a). To order (k F a) 2 the result is id. Il3l Il4|] 
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2 + 07-1) (J^(k F a) + ^(11 - 21og(2))(^a) 2 ) + 0((k F a)^ 
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Effective range corrections appear at 0((kpa) 2 (kpr)) and if g is bigger than 2 logarithmic 
terms appear at 0((/cfo) 4 log(fcpa)). 

This expansion is clearly useless if {kpa) ^> 1. In the case of zero baryon density it is 
well known that an infinite set of bubble diagrams with the leading order contact interaction 
must be resummed if the two-body scattering length is large. It is natural to extend this 
calculation to non-zero baryon density by summing particle-particle bubbles with the finite 
density propagator given in equ. (0J). In traditional nuclear physics this approach is known 
as Brueckner theory 0,13]. The elementary particle-particle bubble is given by 



d 3 q 9g p k F 



■fpp(*,s) . (6) 



(2tt) 3 k 2 -q 2 + ie 4tt (2tt) 2 - 

We are following here the notation of Steele The theta function 9^ = 9{ki — k F )9{k2 — k F ) 
with k\ i = P/2±k requires both momenta to be above the Fermi surface. The first term on 
the RHS is the vacuum contribution which contains the PDS renormalization scale p. The 
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FIG. 2: Interaction energy per particle from particle-particle ladder diagrams as a function of 
(fcj?a). The curves show a numerical calculation of the ladder sum and the two approximation 
fpp — > 2 and /pp -> {fpp) discussed in the text. 

second term is the medium contribution which depends on the scaled relative momentum 
K— k/kp and center-of-mass momentum s — P/(2kp). For s < 1 we have 



fp P (K,s) = 1 + s + /clog 
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Particle-particle ladder diagrams built from the elementary loop integral given in equ. © 
form a geometric series. The contribution of ladder diagrams to the energy per particle is 



given by 



21] 



Epp 3(g - 1)tt 2 r d 3 P d 3 k 
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(8) 



This result can be interpreted as the trace of the in-medium particle-particle scattering 
matrix over all occupied (hole) states. Note that equ. (JBJ) is independent of the renormal- 
ization scale parameter /i. This is in contrast to the perturbative result equ. (0) which is 
independent of fi only if a is small. In general the integral in equ. (jHJ) has to be performed 
numerically. Steele suggested that in the large D limit the function fpp can be replaced 
by its asymptotic value 2 and we will examine this claim in Section Another possible 
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approximation is to replace / PP by its phase space average 

(f PP ) = A(n_ 2 log(2)), (9) 

where (.) denotes an average over all momenta corresponding to occupied (hole) states. In 
this case we find 

£pp_, 2(fc P a)/(3vr) 
A 15 J 2Ml-3| F (ll-21og(2))(^a)- 1 ) 

This approximation has the virtue that the energy per particle from ladder diagrams agrees 

with the perturbative result up to 0((kpa) 2 ). In Fig. |2]we compare the two approximations 

with numerical results. We observe that all calculations agree fairly well if the scattering 

length is either negative or positive and large. For g = 2 the parameter £ in equ. (JIJ is given 

by 

£ = 0.44 (fpp -2), £ = 0.32 (/ PP - (/ PP )), £ = 0.24 (num). (11) 

For g > 2 the results indicate that £ is negative and the homogeneous low density phase 
is unstable. The different calculations shown in Fig. |21 disagree strongly in the regime 
(kpa) ~ 1. Approximating / PP by a constant leads to a singularity in the energy per 
particle. In the numerical calculation this singularity is smoothed out, but a significant 
enhancement in the energy per particle remains. However, even in this case the particle- 
particle ladder sum has singularities for certain momenta that correspond to occupied states. 
These singularities are presumably related to the existence of deeply bound two-body states 
in the vacuum for a ~ In this case interactions between the bound states are essential 
and the approximations used in this section are not reliable. 

III. EFFECTIVE RANGE CORRECTIONS AND HOLE LADDERS 

In this section we study the question whether the ladder sum can be systematically 
improved by including higher order terms in the effective lagrangian. The two-derivative 
term proportional to Ci incorporates effective range corrections. At zero baryon density the 
particle-particle scattering amplitude is [^(J 

T(k) - C ° + ^ (12) 

The effective range approximation corresponds to keeping Co to all orders but treating 
C 2 as a perturbation. The structure of equ. (fT^j) is very simple because in dimensional 
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FIG. 3: Equation of state of a dilute Fermi gas in the unitary limit (fcj?a) -tooasa function of 
the effective range. We show the parameter £ defined in the text as a function of (kpr) for two 
different values of the Fermi momentum in units of the PDS renormalization scale fi. 

regularization, in both the MS and PDS renormalization schemes, powers of momentum 
internal to the elementary particle-particle loop diagrams are converted into powers of the 
external momentum. At finite density the integrals are more complicated. We discuss the 
calculation of the bubble sum in the appendix and only present the results here. The energy 
per particle is 
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where /„ and g n are 




(14) 



(15) 



and satisfy the following relations 



I (k F ,K,s) 



(16) 



I 2 (k F ,n,s) = k 2 I (k F ,K,s) + g 2 (k F ,s), 
I 4 (k F ,K,s) = k 2 I 2 (k F , k, s) + g 4 (k F , s). 



(17) 
(18) 



The explicit forms of g 2 and g± are 
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(19) 
(20) 



and both g 2 and #4 vanish as k F — > 0. If fcpC2 <C Co then equ. (flMj) can be simplified to 



£ PP 3(0 - l)vr 2 r d 3 P d 3 k 



A kp J (2tt) 3 (2tt) 3 ^ fe 1 - C MI 

Using equ. (jHJ) to relate the coupling constants Co and C 2 to the scattering length and the 
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effective range we find 
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We observe that the energy per particle depends on the renormalization scale fi. We are 
particularly interested in the situation when \k F a\ ^> 1 and \k F r\ < 1 for which 
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(23) 



We find that effective range corrections are small and independent of /i provided \k F r\ < 1 
and k F / fi < 1. Evaluating the integral by replacing /pp and #2 by their phase space averages 
and taking k F a — > 00 we get 



£(k F r) = 0.32 + 0.19(k F r) + 0(k F /fi, (k F r) 2 ). 



(24) 



Using equ. (j!3)l we can also study the behavior of the universal parameter £ for larger values 
of k F r. The result is shown in Fig. 01 We observe that the dependence of £ on k F r becomes 
weaker as k F r grows. Jn the limit k F r —>■ 00 the parameter £ slowly approaches the free 
Fermi gas value £ = 1 
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FIG. 4: Energy per particle from hole-hole ladder diagrams as a function of Cq = MkpCo/ {Air) 



Cq is equal to {kpa) in the MS scheme, but in the PDS scheme the relation between Cq and (kpa) 
depends on the renormalization scale fi. We also show the perturbative result up to order Cq. 



In addition to higher order corrections to the effective interaction we also consider larger 
classes of diagrams. A simple extension of the calculation of the particle-particle ladder sum 
is the inclusion of hole-hole ladders. Following the steps that lead to equ. (jHJ) we find 



E H h 3(0 - 1)tt 2 r d 3 P d 3 k 
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f HH (K,s) (25) 



where s) = fpp(K,—s) is the hole-hole bubble. We have subtracted the first two 

terms in the expansion of the geometric series. These terms have UV divergences and need 
to be treated separately. In our case this is not necessary since the two contributions are 
already included in the particle-particle ladder sum. We observe that the remaining part 
of the hole ladders is finite and only depends on Cq and not the PDS renormalization scale 
fi. This implies that if the coupling constant is related to the scattering length according 
to equ. Q, the energy per particle will depend on the renormalization scale /i. Numerical 
results for the hole-ladders are shown in Fig. 0] We observe that if Cq is of natural size 
the energy per particle from hole ladders is indeed very small compared to the contribution 
from particle ladders. 
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IV. PARTICLE-HOLE RING DIAGRAMS AND THE LARGE g EXPANSION 



Another important class of diagrams is the set of particle-hole ring diagrams. In gauge 
theories ring diagrams play a crucial role since they incorporate screening corrections and 
their inclusion of is necessary to achieve a well behaved perturbative expansion. In theories 
with short range interactions particle-hole bubbles also provide important corrections to the 
effective interaction. For example, particle-hole screening corrections reduce the s-wave BCS 
gap by a factor ~ 1/2 in the weak coupling limit. 



The real and imaginary parts of the particle-hole bubble are given by 
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(26) 



(27) 



< v < q - 2- 



q < 2, - . H 2 , 

where v = k§M/k 2 F and q — \k\/kp. Ring diagrams containing particle- hole bubbles can be 
summed in essentially the same way as the particle-particle ladders. The main difference 
arises from different spin and symmetry factors. The spin factor of the n-th order particle- 
particle ladder contribution is g(g — 1)2™. The spin factors of particle- hole diagrams are 
more complicated. The situation simplifies in the limit of large g, often called the large N 
limit, as it is equivalent to the limit of a large number N of degenerate spin 1/2 fermions. 
In this case every particle-hole bubble contributes a factor g. Indeed, one can show that the 
particle- hole ring diagrams are the leading diagrams in the large g limit jlo| . 

The n-th order diagrams in both the particle-particle ladder sum and the particle-hole 
ring sum have symmetry factors 1/n. In the case of the ladder diagrams this factor is 
canceled by a factor n specifying the n different ways in which the diagram can be cut to 
represent it as a particle-particle Green function integrated over all occupied states. For the 
ring diagrams it is more convenient to carry out the energy integration explicitly, and no 
factor n appears. As a consequence, the sum of all ring diagrams is a logarithm. In the 
large g limit we find 
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FIG. 5: Energy per particle from particle-hole ring diagrams as a function of Co = MUfCq/ '(4.7r) 



Co is equal to (fcpci) in the MS scheme, but in the PDS scheme the relation between Co and (kpa) 
depends on the renormalization scale ft. We also show the perturbative result up to order Cq. 

where we have subtracted the first two terms in the expansion of the logarithm in order to 
make the integral convergent. These two terms can be computed separately and correspond 
to the first two terms of the perturbative expansion given in equ. JHJ). Equ. (|28|) shows that 
the correct way to take the large g limit is to keep gCo constant as g — > oo. In this case 
the free Fermi gas contribution as well as the O(kpa) correction in equ. (0) is of order 0(1) 
while the ring diagrams give a correction of order 0(1/ g). 



Since the energy is real, the integral in equ. (J2SJ) can be written as 

E _3 
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SrCoImIIo(i/,g) + (gC ) Imn (^, q)ReU (u, q) 
5f(7oImno(z/,g) N 
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In the PDS scheme Cq is related to the scattering length by equ. (J3J) and the ring energy 
depends of the renormalization scale. Numerical results (for g = 2) are shown in Fig. |3J 
For simplicity we have taken fx = 0. In this case the energy per particle goes to infinity as 
(kpa) — > oo. For other values of fi the energy is finite, but strongly dependent on fi. We also 
observe that the ring energy per particle is less than 3E F /5 for (hp a) < 0.9, which implies 
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an instability of the homogeneous system. 

We emphasized above that equ. (|2*Hj) is correct only in the large g limit. We can also 
compute the ring sum for g = 2. In this case there are two possible channels with total spin 
zero and one. The n-th order ring diagram has spin factor (—2, 4, —2, 4, . . .). The sum of all 
ring diagrams is 

E = - l -l-^[log(l-C n(k))+3log(l + C U(k)) 



-2C U (k)+2(C U (k)) 2 . (30) 



The integral can be calculated as in equ. (J29|) and the result is shown in Fig. We observe 
that although the correct g = 2 result is quite different from the g — > oo result evaluated 
at g = 2, qualitative features, like E/A — > oo as Co - > oo and the presence of an unstable 
regime with Ej A < 3Ep/5, remain unchanged. 



V. LARGE D EXPANSION 



In the previous section we saw that the particle- hole ring energy can be interpreted as the 
leading order contribution to the energy in the large g limit. This raises the question whether 
there exists an expansion that gives the ladder sum as the leading order contribution. In 
a very interesting paper Steele suggested that expansing in 1/D, where D is the number 
of space-time dimensions, is the desired scheme Q. If true the 1/D expansion offers a 
systematic approach to the fermion many-body system in the limit (kpa) — > oo. 

The main idea is that many body diagrams in a degenerate Fermi system are very sensitive 
to the available phase space and that the scaling behavior of phase space factors in the large 
D limit could be a basis for a geometric expansion. Steele argued that the 1/D expansion 
corresponds to the hole-line expansion in traditional nuclear physics, and that it is consistent 
with EFT power counting for systems with a large scattering length. In this section we shall 
examine these claims in more detail. 

The spatial density of a free Fermi gas in D space-time dimensions is given by 

P ~ (2n)^D-V " D - 1 ~ r(^i) ' (31) 
where is the surface area of a D-dimensional unit ball. In the following we will always 



12 




0-0D= 


=4 


B-H D= 


=10 


<^>D= 


=20 




=30 


<3-<] D= 


=40 




0.3 0.4 
s,k [kp] 



FIG. 6: Scaled particle-particle scattering amplitude Dfpp(K,s) for k = s (left panel) and s = 
(right panel). We observe that except in the BCS limit s = 0, k — > 1 the function Dfpp approaches 
a smooth limit as D — > oo. 



take the degeneracy factor g = 2. The energy per particle is 
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We can compute perturbative corrections to this result in D dimensions. To leading order 
in Co we find 
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This expression indicates that in the weak coupling limit the large D limit should be taken 
in such a way that 



A = 



n D - l c kp'- s M 



const. 



(34) 



D(2n) D - 1 

In the following we wish to study whether this limit is smooth even if the theory is non- 
perturbative, and what class of diagrams is dominant. For this purpose we consider the in 
medium particle-particle bubble for an arbitrary number of space-time dimensions D. The 
result is 

d D - x q 



(27r) D - x k 2 - q 2 + ie 



fvac{k) + 



where 
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with xo(s,t) = (s 2 + t 2 — \)/(2st). The factor 2 was inserted in equ. (j^Bj) so that the 
normalization of fpp is consistent with the previous case D = 4. The other terms are 

c(D, x ) = 2x 2 F l Q, 2 - y , | , (37) 
c D = c(D,l) = V g -n ' (38) 



2^1 (l, 



2 

D-l D+l !- 



= p v _ , 1) ^( 1 +; ) i-^ > ( 39 ) 

where 2-^1 (a, 6, c, 2) is the hypergeometric function. Numerical results for the function 
fpp(s, k) are shown in Fig. El We observe that the particle-particle bubble scales as 

f P D p\s, k) = ^f PP (s, K)(l + (i)) . (40) 
The scaling behavior can be verified analytically in certain limits. We find, in particular, 

/ao.o)—. (41) 

There is a subtlety associated with the BCS singularity at s = 0, k — > 1. Fig. shows that 
the logarithmic singularity is not suppressed by We observe, however, that the range 

of momenta for which fpp is enhanced shrinks to zero as D — » 00. We will study pairing in 
the large Z) limit in Sect. IVI1 

Since f p D p ~ 1/Z? we conclude that if the large D limit is taken according to equ. (|3~3jl 
then all ladder diagrams with particle-particle bubbles are of the same order in 1/D. The 
sum of all ladder diagrams can be calculated by noting that, except for the logarithmic 
(BCS) singularity at s = 0, k = 1, the particle-particle bubble is a smooth function of the 
kinematic variables s and k. Hole-hole phase space, on the other hand, is strongly peaked 
at s = R = l/y/2 in the large D limit. This can be seen by re-expressing s, k in terms of 
and using A; 12 — > kp as D —> 00. If the phase space is strongly peaked we can replace the 
function fpp(s, k) by its value at s, R. We have not been able to calculate the large D limit 
of fpp(s,R) analytically. Our numerical results show that limD^oo fpp (s,R)/ fpp (0,0) = 
2.02 ± 0.02. We therefore conjecture that /^(s, R) = 1/D ■ (1 + 0(l/D))). 

We illustrate the method by calculating the second order correction. This contribution 
involves an integral of the particle-particle bubble over hole-hole phase space. We find 



dP^P f d^k p) ifcjP- 2 
(2n) D -iJ ( 2 7r) D -^ kIpp[K ' S) - D 2 
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and the energy per particle is given by 

n D -iC k$- 3 M] 2 [ k 2 F \ (43) 



i = 2 
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Higher order terms can be calculated in the same fashion. As in D = 4, particle-particle 
ladder diagrams sum to a geometric series. We get 

E L A „ / 1 \1 / k 2 F ^ 



a r 1-2A {d)]{w)' < 44 > 

where A is the coupling constant defined in equ. (jHU). We observe that if the strong coupling 
limit A — > oo is taken after the limit D — > oo the universal parameter £ is given by 1/2. 

We have been able to compute the particle ladder contribution to the energy per particle 
in the large D limit. Steele argued that all other contributions are suppressed by powers of 
1/D since each additional hole line involves an integral over the Fermi surface of the type 
shown in equ. ()31|) which gives a least one power of 1/D. It is not clear if this argument 
is entirely correct. We have found, for example, that the main contribution of the particle- 
particle bubble also scales as 1/D. We study this problem in more detail in Sees. IVIIand lVIII 
where the potentially relevant pairing and screening corrections are examined, respectively. 

VI. PAIRING IN THE LARGE D LIMIT 

In the previous section we noticed that in the large D limit the particle-particle bubble is 
enhanced when s = and k — > 1. In this limit the two particles are on opposite sides of the 
Fermi surface and the logarithmic enhancement of fpp is the well known BCS singularity. 
The result suggests that in the large D limit the pairing energy might dominate all other 
contributions to the energy density. In this section we shall study this question by computing 
the BCS gap and the pairing energy in the large D limit. 

If the interaction is weak and attractive we can derive the standard BCS gap equation 



(see, for example, reference 



m 



\C \ r d D ^p A 



with e p = E p - Ep and E^ = »V(2M). The integral in equ. ()45|) can be carried analytically 
for arbitrary D. We find |25l l26j 



DXty 2 W 2 



sin(7ra) ^ 1 \ a/I + x 



1 = ^7rrr(i + ^j M-TTTT^) (J(,) 
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where A is the dimensionless coupling constant defined in equ. (|34|i. x = A/Ep is the 
dimensionless gap, P a (z) is the Legendre function and a = (D — 3)/2. We note that by 
going to arbitrary D we have regularized the UV divergence in the gap equation using 
dimensional regularization. If the gap is small, x< 1, equ. (J46j) can be solved using the 
asymptotic behavior of the Legendre function P a (z) near the logarithmic singularity at 

sinfavr) f, (1 + z s 



P a {z) ^ (log f — — ) + 27 + 2^(a + 1) + 7TCOt(«7r)J . 



(47) 



To leading order in 1/D we can also use the asymptotic expression for the Digamma function 
ip(a) = T'(a)/T(a) ~ log(a) + 0(l/a). We find 



exp 



(48) 



D V D\, 

where 7 ~ 0.5772 is Euler's constant. We observe that the exponential suppression of the 
gap disappears if the large D limit is taken at fixed A. However, the exponential suppression 
in A is replaced by a power suppression in 1/D. 



Next we ca 
limit we have 



culate the pairing contribution to the energy density. In the weak coupling 
2a 



E 



d 



D-l 



P 



A 2 



(2< 



D-l 



2je 2 ; + A2 



+ J e^ p + A 2 - e p 



(49) 



The integrals can be calculated in the same fashion as the integral that appears in the gap 
equation. We find 



E 



D-l 



(2< 



D- 



r ^"^(l + x 2 ) 



a/2 



1 



ax 



a: 



x- 



P, 



2 2a + 4 
/ 1 



1 



X* 



a+1 



(50) 



a + 2 UT1 v VTT^ 
In the limit x — >• the logarithmic singularities in the two terms in the curly brackets cancel 
and the energy per particle is proportional to x 2 . We find 

D 



E 
A 



— Ef \e;. 



(51) 



Since A/Ep = 0(1/ D) we conclude that the pairing energy per particle scales as 1/D 
in the large D limit. This implies that the pairing energy is suppressed compared to the 
contribution from the ladder sum given in equ. 
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FIG. 7: First order screening correction to the effective particle-particle interaction. 



VII. SCREENING IN THE LARGE D LIMIT 

In D = 4 dimensions screening of the elementary four fermion interaction by particle- 
hole pairs gives an important contribution to the effective interaction. It is well known, 
for example, that screening reduces the magnitude of the gap in the weak coupling limit 
by a factor (4e) 1//3 ~ 2.2 j29|. On the other hand, if the 1/D expansion corresponds to an 
expansion in the number of hole lines then we expect that the screening correction should 
scale as 1/D in the large D limit. The basic particle-hole bubble is given by 

n(f, q) = Tg^pzi fPH\v, q) (52) 

with 

f PB {u, q) = - t k D ~ 2 dk f 1 dx(l - xY /2 - 2 ^^i (53) 
c D Jo J-i v l - u z qk 

where u q k = qkx + q 2 /2. This integral can be evaluated analytically in D = 4 and reduces 

to equ. (|26j). The particle- hole bubble leads to a renormalization of the effective interaction 

as shown in Fig. [7| In the weak coupling limit only the interaction of two quasi-particles 

near the Fermi surface is important. For s-wave pairing we can write 

where fpu is an average over the Fermi surface 

Iph = ^ J\ dx (1 - x 2 ) D/2 ' 2 fpH (o, y/2{l - x)) • (55) 

Replacing the bare interaction with the effective one in the BCS gap equation leads to a 

correction for the pairing gap. We find 

2e~^E F / — \ / 1 \ 
A= — - — exp (-/p H J exp (-— j . (56) 

For D = 4 the integral in equ. (|55|l can also be evaluated analytically, yielding f PH = 
(21og(2) + l)/3 ~ 0.79. Numerical results for D > 4 are shown in Fig. |H| We observe that 
the screening correction vanishes as 1/D for large D. 
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FIG. 8: Screening correction to the effective s-wave particle-particle interaction as a function of 
the number of space-time dimensions D. 

VIII. CONCLUSIONS 

In this work we considered different many body theories for a system of non-relativistic 
fermions described by an effective field theory. We are interested, in particular, in a sys- 
tematic approach to the problem of a dilute liquid of fermions in the limit in which the 
scattering length is large compared to the inter-particle spacing. This led us to consider the 
large g and large D expansions. 

At leading order the large g expansion corresponds to summing all particle-hole ring 
diagrams with the additional approximation that all spin factors are replaced by their large g 
limits. The problem is that the large g limit is not suitable for studying the limit \kpa\ — > oo. 
Indeed, the naive large g limit corresponds to taking \kpa\ — > 0. This is also manifest in 
the functional form of E/A at leading or sub-leading order in 1/g. The energy per particle 
in the limit \kpa\ — > oo is strongly dependent on the regularization scale. The only sensible 
alternative is the Bose limit where we keep kpa fixed and small and take g — > oo with 
p = gk F /(6n 2 ) constant 

From the study of the two-body system in effective field theory we know that if the 
scattering length is large then the two-body interaction has to summed to all orders. This 
suggests that a sensible many body theory has to contain at least all particle ladders. We 



18 



E9 





i 1 


1 'ill 1 

O i 




o pp ladders 


i cj 




- pp(f pp ->2) 






-— PP(fpp->fp P ) 


; o, _ 

oi 1 




% pp (D -> oo) 


o t 

* 0\ 






\o* 
o \o\ 

*0 — 


BCS 


o 




crossover 






o 

X* o 

K\ o° 
\ o 


1 , 1 ', >, 1 



l/k„a 

F 



FIG. 9: Total energy of an interacting fermion gas in units of the energy of a free fermion gas as a 
function of (A^a) -1 . The open circles show the result of a numerical calculation of particle-particle 
ladder diagrams. The dashed and dash-dotted curves show the two approximations discussed in 
Sect, [n] The star is the result of the D — > oo calculation in the unitary limit. 



have studied the particle ladder sum in effective field theory. If the particle-particle bubble 
is replaced by its phase space average then the particle sum can be carried out analytically. 
This approximation gives £ = 0.32 for the universal parameter in the equation of state |3Jj. 
Numerically, we find a smaller value £ ~ 0.24. The ladder approximation is unreliable in the 
regime (kpa) ~ 1 in which deeply bound two-body state are important. These results are 
summarized in Fig. El We observe that different approximations agree in the weak coupling 
(BCS) limit, but give a range of predictions for the crossover behavior. 

It is clearly desirable to construct a systematic expansion that contains the ladder sum 
at leading order. In a very interesting paper Steele suggested that the large D expansion 
provides the desired approximation scheme. In order to test this idea we have calculated 
particle ladders for arbitrary D. We find that if the coupling constant is scaled appropriately 
then all particle ladders are indeed of the same order in 1/D. We also find that the universal 
parameter £ is given by £ = 1/2, in surprisingly good agreement with the Green function 
Monte Carlo result £ = 0.44 22|. We also verified that the contribution from the pairing 
energy is suppressed by 1/D. 

There are many questions regarding the 1/D expansion that remain to be addressed. We 
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have not been able to construct a general method for computing 1/D corrections. We have 
also not succeeded in showing that the 1/D expansion corresponds to an expansion in the 
number of hole lines. Our result for the energy per particle in the large D limit differs from 
Steele's result £ = 4/9 because he did not compute many body diagrams for a general number 
of space-time dimensions. Instead, he performed certain kinematic expansions in the D = 4 
loop integrals. Even if this method was correct it would correspond to a partial resummation 
of 1/D corrections. We also believe that the expansions that are used in Steele's paper are 
not convergent. Finally, we note that there is an argument due to Nussinov and Nussinov 
[3] which indicates that for D > 4 the ground state consists of non-interacting, zero energy 
bosons and that £ = 0. This argument may indicate that there is a subtlety with regard to 
the order of the D — ► oo and (kpa) — > oo limits. 
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APPENDIX A: EFFECTIVE RANGE CORRECTIONS 

In this appendix we provide some details regarding the calculation of effective range 
corrections to the particle-particle ladder sum. The effective interaction is given by 

(P/2 ± h\V\P/2 ± k 2 ) = C + ^-(kl + k 2 2 ) = C + C 2 L • k\ + C* • k 2 2 , (Al) 

where C% = = We consider the following amplitudes 

• X§(X^): The sum of all amplitudes T(P/2 ± k, P/2 ±fc) starting with a R(L) vertex 
and containing N Gi (L or R) vertices. 

• X^'- The sum of all amplitudes starting with a Co vertex and containing N C2 (L or 
R) vertices. 
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FIG. 10: The recursion relations between A N ,A N and A N . 

In this appendix we evaluate Y^n=T(^n + X§ + -^at)- O ur strategy is to derive a set of 
recursion relations between the amplitudes and then use these relations to compute the sum. 
We define 

-R _ aR v0 _ ao (A2) 



Y L — h 2 ■ A L X u — A 1 



X° - A 



The only difference between Xn and An is that An does not include the contribution from 
the initial pair momentum k. The amplitudes An satisfy (see Fig. irU|) 



A L - C<2 T A M C2 T A R M C2 T A° 
A N+i - -y-^A/V.L + —loA N + —1qA n , 



A R C<2 T A -u ^ T A R _i_ ^ T ,1" 

A N+l - -^UAn,l + —i 2 A N + —1 2 A N , 

" C 



(A3) 



A% = V -^ T [hA L N + I A N ]. 

Since A° N can be expressed as the combination of A N and A N one can simplify the above 
recursion relations and obtain 



1 - C Io)A N+1 



(1 - C I )A N+1 



L 


C 2 r 


N+l — 


2 - 


R 


C 2 r 


N+l — 


2 - 



I 2 A L N + I A N 



(h - C Ioh + C li)A L N + I 2 A 



(A4) 



The crucial point is that the Ihs of equ. (|A4J) starts from A 2 while the rhs starts from A\. 
Therefore one obtains 

/ oo \ oo oo 

2(1 -Colo) [T, A* -Af) = C 2 I 2 £ A L N + C 2 I £ A N , 



\N=1 



N=l 
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N=l 



2(1 -Colo) 

\w=i 

where and Af are given by 
A? 



R 



oo oo 

C 2 (h - C I h + C Q ll) Y. A n + C 2 h E A m (A5) 

N=l N=l 



C 2 C2 Co _ C2 

T + T ' 1 - C / ~ 2(l-C / )' 



— £- 2 

2 fc 



Co 



Co 



'2 • 



2" 1-Co/o 

We define the sums over ^4^r, ^v ^^d A^ as 

00 00 

J2 a n = y r , E ^ 



2 



/c 2 + C (I 2 - k 2 I ) 

1 - c / 



(A6) 



N=l 



N=l 



n, x: a °n = y o 



(A7) 



AT=l 



After some algebra one obtains 

C 2 (l-C / )(l-C 2 (/ 2 -A: 2 /o)/2) 



L ~ ~2 ' (1 - C /o - C 2 / 2 /2)2 - C 2 / (Co(-/o/4 + if) + / 4 )/4' 

. C2 (1 - Cq/q - C 2 / 2 /2)(fc 2 + C (/ 2 - fc 2 / )) + c 2 (c (-y + g) + / 4 )/2 
o ' n _ n. r_ _ r_ /o^2 nijjnj r_ r. 1 r2N , r.W/i d A °J 



(1 - Cq/q - C 2 / 2 /2) 2 - C|/o(Co(-/ /4 + If) + I 4 )/4 



It is now straightforward to evaluate the sums over Xjv 



N=oo 



E (X^ + X* + X° N ) = k 2 - Y L + Y R + Y 

N=l 



Cn/. 



± Yr + ( k 2 + 02 

Co/o V 1 — Co/| 

Co (1 — Cn/n — Co- 



(1 - Cq/q - C 2 / 2 /2)(fc 2 + C (/ 2 - k 2 I )) + C 2 (C (- V 4 + I 2 ) + J 4 )/2 

2(1 - Coh) ' (1 - C /o - C 2 / 2 /2) 2 - C 2 2 J (Co(-/o/4 + I 2 2 ) + / 4 )/4 

C 2 ((A: 2 + C (/ 2 - fc 2 / ))(l - C 2 (/ 2 - fc 2 J )/2) 



~2~" (1- C /o - C 2 / 2 /2)2 - C 2 / (Co(-/o/4 + if) + / 4 )/4 
Equ. (|A9J) can be used to derive equ. (fT3j) . 
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